Phase diagram of magnetic configurations for soft 
magnetic nanodots of circular and elliptical shape 
obtained by micromagnetic simulation 

E.R.P. Novai£] and A. P. Guimarae^ 

Centra Brasileiro de Pesquisas Fisicas, 
R. Xavier Sigaud 150, 22290-180, Rio de Janeiro, Brazil. 

Abstract 

Magnetic disks or dots of soft magnetic material of sub-micron dimensions may have as the lowest 
energy magnetic configuration a single-domain structure, with magnetization either perpendicular of 
parallel to the plane, or else may form magnetic vortices. The properties of these vortices may be used 
to encode data bits, in magnetic memory applications. In the present work the OOMMF code was used 
to compute by micromagnetic simulation the energy and the magnetization of circular and elliptical 
nanodots of permalloy. For the elliptical magnetic dots the analysis was made for variable thickness and 
length of major axis, keeping a 2:1 axis ratio. From the simulations, a phase diagram was constructed, 
where the ground state configurations of the nanodots are represented in a diagram of nanodot height 
versus length of the major axis 2a of the ellipse. The phase diagram obtained includes regions with one 
and two vortices; it is similar, but more complex than that derived using a numerical scaling approach, 
since it includes configurations with lateral vortices. These diagrams are useful as guides for the choice 
of dimensions of elliptical nanodots for practical applications. 

1 Introduction 

Nanoscopic and mcsoscopic magnetic structures in the form of dots of circular or elliptical shape have 
attracted the interest of many workers in recent years in view of their very interesting physical properties 
and for their potential applications; for an introduction see [1]. Quasi-twodimensional magnetic objects 
of nanoscopic or sub-micron dimensions (referred here as nanodots) made of soft magnetic materials, such 
as permalloy, may present, for their lowest energy configuration, a single-domain ordering, with magnetic 
moments either perpendicular or parallel to the plane of the dot. For single-domain, we understand the 
magnetic structure where the spins are in general parallel, although deviations of the spins near the surfaces 
are expected. They may also, depending on the dimensions, form magnetic vortices, or swirls. 

These magnetic vortices appear as lowest energy magnetization configurations and have been observed 
by many experimental techniques, such as magnetic force microscopy [5], X-ray microscopy [5], or inferred 
from hysteresis curves |17| . The magnetic vortices have a core that is magnetized perpendicularly to the 
plane of the nanodot, with a radius of the order of the exchange length in the thin dot limit [H, about 5 
nm in the case of permalloy. Applied magnetic fields or spin polarized currents may reverse the direction of 
magnetization of this core; otherwise, the magnetic vortices are stable structures. 

The proposed applications of magnetic nanodots include their use in patterned magnetic recording media 
[5] and as elements in magneto-resistance random access memories (MRAM's) [6], [7|. Nanodots of circular 
shape present flux closure, and this implies that the perturbation from neighbor elements is reduced, allowing 
the use of higher density arrays, an important consideration in applications. 

Nanodots exhibit different magnetic structures, depending on their dimensions and geometry. Consider, 
for example, 20 nm-thick circular nanodots prepared from a soft magnetic material, such as permalloy. Below 
a radius of 10 nm, the magnetization is perpendicular to the plane, from 10 to about 30 nm the magnetization 
is parallel to the plane, and above this radius, the minimum energy magnetic structure is a vortex. 
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Magnetic vortices are structures characterized by a direction of rotation (circulation c), that may be 
clockwise (CW, c = —1) or counterclockwise (CCW, c = +1), the direction of the magnetization of the vortex 
core, that may be along the positive z axis or antiparallel to it; this defines the polarity as p = +1 or p = —1, 
respectively. The core has a radius of the order of the exchange length of the material l^^ = ■\/2y4//ioA/| 
(or = ^2j4/47rA/| in the CGS), where A is the exchange stiffness constant and Ms is the saturation 
magnetization. With the parameters used for permalloy in the present work, = 5.29 nm. The combination 
of circulation and polarity defines the chirality, or handedness of the vortex; a vortex may be righthanded 
[cp = +1) or Icfthandcd (cp = — 1); one should note that the word chirality is sometimes used in the literature 
to describe the sense of rotation of the vortex. 

The different magnetic configurations, and consequently the large variation in the magnetic properties 
observed in nanodots as a function of dimensions, underlines the interest in the study of diagrams (phase 
diagrams) mapping the parameter space where a given magnetic behavior is to be expected. We have chosen 
in the present work to draw such diagrams for nanodots of circular and elliptical shape using micromagnetic 
simulation. This simulation allows the computation of the total energy of the nanodot, which encompasses a 
magnctostatic contribution, the exchange contribution and the anisotropy term. Micromagnetic simulations 
are invaluable for the prediction of magnetic behavior of these nanostructures, although necessarily containing 
some idealizations, for example, assuming the constancy of the thickness of the dots, homogenous physical 
properties, or the absence of defects. 

In the present work the energy per unit volume and the equilibrium magnetic structures of thin permalloy 
islands of circular or elliptical shape and sub-micron dimensions were studied using the micromagnetic 
simulation program OOMMF, a free software available from NIST 

citeoommf. The crystalline anisotropy term in the energy was neglected, as it is usual in simulations of 
permalloy samples. 

2 Results and discussion 

We have made a systematic study of the magnetic configuration of circular and elliptical (of axial ratio 2:1) 
nanodots of permalloy through micromagnetic simulations. For the simulations we used the 3D OOMMF 
code. The computation was performed for permalloy samples, using the exchange stiffness constant A = 
1.3 X 10~^^ J/m, the saturation magnetization Mg = 860 x 10'^ A/m and the anisotropy constant if = 0, the 
typical values used in the program. To make the present results of more general use, they have been given 
in terms of parameters normalized using the exchange length l^x- The cell size in most cases was taken as 
5x5x5 nm'^; the dependence of the simulation results on the cell size was verified for some configurations. 

For some dimensions of the nanodots. the simulation may converge to a configuration that docs not 
correspond to the absolute energy minimum. Therefore, the simulations made with parameters near the 
boundary regions between different configurations had to be made by first imposing every possible config- 
uration, and afterwards comparing the resulting energies to determine the spin arrangement corresponding 
to the absolute energy minimum. 

To save computation time, the simulation for the larger elliptical dots used as a starting configuration a 
multidomain structure that was close to a single or double vortex structure, depending on the value of the 
major axis of the ellipse whose magnetic structure was being simulated. Some checks with the full calculation 
were made to verify that this procedure did not lead to secondary energy minima. 

2.1 Circular nanodots 

Based on the shape of the hysteresis curves of nanodisks of Ni, Co, CoNi and CoP of different thicknesses 
and diameters, Ross et al. [8] have drawn a magnetic configuration phase diagram. The diagram con- 
tained three configurations, described as in-plane flower, out of plane flower, and vortex/multidomain. The 
boundaries between the regions corresponding to these arrangements agreed with the frontiers computed by 
micromagnetic simulation. 



2 



Metlov and Guslienko [5] derived an expression for the total energy of circular nanodots, and from this 
result obtained a phase diagram with regions of in-plane magnetization, perpendicular magnetization and 
vortex structure. Scholz et al. |10] computed the equilibrium magnetic configuration for some dimensions of 
dots by micromagnetic simulation and found a general agreement with the phase diagram derived by Metlov 
and Guslienko [S]. 

Another simulation study, this time using a scaling approach, was made for circular and elliptical nanodots 
[llj . For circular dots the phase diagram presented in that work contained parallel, perpendicular and vortex 
magnetization structures; their results for elliptical dots will be commented below. 

Our results for the total energy per unit volume {E/V), or energy density, for some dimensions of circular 
permalloy nanodots are presented in Figs. [T]and[ll The first figure shows the variation of E/V versus reduced 
radius r/l^x for a 20 nm thick nanodot, for three magnetization configurations; it is evident that the energy 
of the minimum energy configuration changes only by a factor of about 3 in the range of radii shown. The 
same is valid for the second figure, where E/V is drawn against reduced thickness h/l^x- The range of energy 
per unit volume for each given magnetic configuration is wider, up to about 10 times in some cases, for the 
dimensions shown in the figures. 




Figure 1: (Color online) Energy per volume versus reduced radius (r/lcx) for 20 nm-thick circular nanodot 
of permalloy, computed by micromagnetic simulation. The quantity lex is the exchange length. 

The computed curves of energy density versus thickness, or height, arc shown, for circular nanodots, in 
Fig. [2^; the height is measured in units of the exchange length l^x- 

Fig. [3] shows the phase diagram for circular magnetic nanodots, obtained from the computed energies 
for the different magnetization configurations. The diagram shows three regions, corresponding to in-plane 
single domain, out of plane single domain, and magnetic vortex. The boundary lines between the several 
magnetic configurations exhibited in Fig. [3] were obtained from the intersection of the curves of E/V versus 
dimension for the different arrangements. 
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Figure 2: (Color online) Energy per unit volume for circular and elliptical permalloy nanodots versus reduced 
height /i//ex for a) 70 nm-diameter circular nanodot and b) elliptical nanodot of major axis 2a = 70 nm. 
The curves were computed with micromagnetic simulation. 



2.2 Elliptical nanodots 

Experimental studies using an array of micrometric elliptical permalloy nanodots with 2:1 axis ratio showed 
both single vortex and double vortex states, the latter configuration being more probable when observed (at 
remanence) after saturation in the direction of the long axis |12] . 

Magnetic configurations with one, two and three vortices were observed experimentally with the MFM 
technique on permalloy ellipses some microns long [13] . In the same study these magnetization arrangements 
were also derived from micromagnetic simulations performed with permalloy ellipses of 720x240x30 nm'^ . 
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Figure 3: (Color online) Phase diagram for circular nanodots of permalloy, as a function of reduced height 
h/lcx and reduced radius r/lcx, drawn from the minimum energy computed with the simulation. The 
continuous line was drawn following a procedure described in the text. 
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Figure 4: (Color online) Phase diagram for elliptical nanodot of permalloy, as a function of the reduced 
height h/lcx and reduced length of the major axis of the ellipses 2a/lcx- The ellipses in every case have axes 
in the ratio 2:1. Note the break in the horizontal scale, to reveal the appearance of two-vortex structures 
obtained for larger ellipses. The continuous line was drawn following a procedure described in the text. 
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Earlier work on elliptical dots using micromagnetic simulation sketched a phase diagram containing three 
configurations: a region of quasi-single domain structure, single domain and vortex. The computation was 
made for fixed nanodot thickness (30 nm) [H] , 

The phase diagram obtained by Zhang et al. [TT] for elliptical dots shows, besides the configurations 
observed for circular dots, a double vortex magnetic arrangement, which appears, for the parameters chosen 
in their illustration (Fig. 7, in for major axes 2a above about 300 nm. In their simulation, the 

contribution of the vortex core to the energy was neglected; the authors estimated that the inclusion of the 
core contribution would displace the boundaries in the phase diagram by about 35 %, quite a significant 
shift. 

The computed curves of energy density versus thickness, or height, are shown, for elliptical nanodots, in 
Fig. [2b; the height is measured in units of the exchange length Ig^- The comparison between the graph for 
circular and elliptical dots shows that the configuration of lateral vortex is a minimum energy configuration 
for elliptical dots, but not for circular dots. In the work of Ha et al. [15] using another computer code, in 
one single case a lateral vortex was observed in simulations for circular dots; however, their methodology 
was different, involving the determination of the magnetization configuration at remanence, after saturating 
the nanodot with an external magnetic field. This difference may explain also other differences in the phase 
diagram that they have obtained. 

The phase diagram obtained from our micromagnetic simulations for elliptical nanodots is shown in Fig. 
m this diagram shows single-domain arrangements, both in-plane and perpendicular, as found in simulations 
for the circular dots. However, it has a more complex structure, with configurations with lateral vortices 
and, for the larger dots, configurations with more than one vortex; the region where two vortices appear is 
shown in the figure. 

We have also shown, in the graph for elliptical nanodots, the curves of energy density for configurations 
of one and two magnetic vortices (Fig. [5|), as a function of dot thickness. One can see from these curves that 
the configuration with more than one vortex, for ellipses with 2a = 2560 nm, corresponds to a minimum in 
energy only for thickness above about 50 nm. 
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Figure 5: (Color online) Energy per volume (E/V) versus thickness h of elliptical nanodot. with axis 2a = 
2560 nm, showing that above about h = 50 nm, the configuration with two vortices has a smaller energy 
density than the configuration with one vortex, from the micromagnetic simulation. 
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Finally, the effects of discretization, inherent in the methodology used here (see [H]), are responsible for 
very small changes in the boundaries of our phase diagram for circular and elliptical nanodots, since on the 
one hand, the reduction of the dimension of the cell for the regions where the vortex corresponds to the 
minimum energy still results in this same configuration as minimum; on the other hand, where the vortex is 
not a minimum this is not altered by choosing a smaller cell size. These effects do not exist for perpendicular 
magnetization. 

3 Conclusions 

Using 3D numerical simulation we have obtained curves of energy per unit volume versus dimension for 
nanodots of elliptical and circular shape that illustrate how energy considerations explain the variability in 
magnetic arrangement found in these magnetic nanostructures. 

We have shown in the present work the regions of geometrical parameters where soft magnetic dots of 
nanoscopic or mesoscopic dimensions, of circular and elliptical shapes, take up single domain and vortex 
magnetization configurations. For circular dots our results are in agreement with previous work performed 
using different methods. For elliptical dots there are some important differences, such as the appearance of 
a region where a magnetic structure with lateral vortices corresponds to an energy minimum. 

Investigations that allow the mapping of these different magnetic configurations are useful in designing ex- 
periments to study the basic properties of these novel magnetic structures, or tailoring them for technological 
applications, such as magnetic random access memories. 
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